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Abstract 



We demonstrate that the neutrino kinetic equation derived by the standard Bo- 
golyubov method is formally gauge non-invariant and give a recipe how to recast it 
to the gauge invariant form recovering the standard Lorentz form weak force term 
which apparently conserves the lepton current. The analogy with the QED plasma 
case is traced which in asymptotic regions implies substitution of electric charge e 
by the induced electric charge e % „ d of neutrino within the phase factor connecting 
CN ■ neutrino gauge invariant and non- invariant distribution functions. 
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^ . 1 Introduction 

«?; 

The neutrino Relativistic Kinetic Equation (RKE) is a useful tool to describe 
many phenomena in astrophysics and cosmology. In particular, neutrinos play 
the most important role for a supernova (SN) burst or in the lepton asymmetry 
formation before the primordial nucleosynthesis in the early universe. The 
usual motivation to use the RKE approach for neutrino propagation in a 
dense matter is stipulated by the account of neutrino collisions: within a SN 
neutrinosphere or in the hot lepton plasma of the early universe before neutrino 
decoupling. 

However, in addition to collision integrals there are self-consistent weak inter- 
action terms in the neutrino RKE [1] that are linear over the Fermi constant 
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~ Gf and analogous to the Lorentz force terms for charge particles in the 
standard Boltzman RKE which in turn are linear over the electric charge ~ q 
(q — — | e | for electrons). 

In the standard kinetics these self-consistent electromagnetic fields are well- 
known to play a very crucial role. In collisionless, or Vlasov approximation, 
such kinetic equations describe, e.g. thermonuclear plasmas in laboratory and 
stars for which an energy exchange between electromagnetic waves (eigen 
modes) and charged particles proceeds faster than via the direct particle col- 
lisions with all following issues in collisionless plasma: instabilities, heating, 
etc. 

One expects that the self-consistent weak interaction (~ G F ) could lead for 
neutrinos to some analogous collective interaction effects, e.g. to neutrino 
driven streaming instability in an isotropic plasma [2], to generation of mag- 
netic fields in the early universe or in a SN [3,4]. 

Recently neutrino RKE has been rederived along different ways in [2,4,5] and 
the goals of this letter are: (i) to establish the conformity of these RKE's 
with the standard Bogolyubov approach (Bogolyubov-Born-Green-Kirkwood- 
Yvon (BBQKY) chains of kinetic equations) used in [1]; (ii) to elucidate the 
physical sense of the ponderomotive weak force appeared in effective Lorentz 
form meanwhile the neutrino density matrix is gauge invariant with respect to 
the corresponding electron current transformation. 

To this end, in section 2 we find that the formal lepton current nonconservation 
in the master RKE [1] stems from the absense of gauge invariance (see eq. (4) 
below) for the usual in quantum statistics definition of the gauge non-invariant 
distribution f(x u x 2 ,t) = Tr(p(f)* + (£ 2 )*(£i)). 

Then in section 3 we suggest the recipe of gauge invariance restoration which 
allows to derive the Lorentz form of the weak force term in RKE [2,5] and 
after that in the main section 4 we find the phase transformation (12) which 
connects the gauge invariant distribution for neutrinos with the gauge non- 
invariant one. 

The above transformation is fully analogous to the one used in the case of 
QED plasma, being based on the same way of inclusion of interaction: with 
the self-consistent e/m field in QED plasma, and with the self-consistent 
electron four-current J M for neutrinos. As a result, these transformations turn 
out to be equivalent with the accuracy of a change of charges e — > e™ d , where 
e™ d is the induced electric charge of neutrino [6] (see section 4) . This analogy 
is completed in section 5 where we recast the weak force term into a form of 
radiative damping force due to a plasmon emission by neutrino. 

In section 6 we give conclusions and in Appendix we remind some well-known 
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properties of self-consistent electromagnetic fields in plasma. 



2 Neutrino RKE for gauge non-invariant Wigner distribution 



Neglecting electron spin from the quantum Liouville equation one finds in the 
Vlasov approximation the neutrino Relativistic Kinetic Equation (RKE) as the 
classical equation for the gauge non-invariant distribution function p v \q, x, t) 

[1], 



±x — + Q ^ = 0, (1) 



dt 



dx 



where x = n = q/q is the velocity of massless neutrino, the derivative q is 
given by 



'q = G F V2c v \Vn^(x,t) - V(nj (e) (f , t)) 



(2) 



Here j( e )(f,t) = (n^(x , t); ] {e \x , t)) = J(d 3 p/(27r) 3 )(p,/e p )f^(p, x, t) 
= (x, t)/e is the four- vector of the electron current density devided on the 
electron charge e = — | e |; cy = 2£ ± 0.5 is the weak vector coupling (upper 
sign for electron neutrinos), £ = sin 2 9 W ~ 0.23 is the Weinberg parameter. 

Obviously, due to the second term in (2) the RKE (1) does not obey the 
neutrino current conservation law, 

djM(x,t) , 

where ffi(x,t) = J d 3 q(q fM / q) (q, x, t) / (2tt) 3 is the neutrino four-current in 
medium. 

First, note that RKE (1) can be derived from the canonical equation dp u ^ /dt+ 
{H, /M} = with use of the neutrino Hamiltonian in a medium [2] 



H = H + V eff = V(Q - G F V2c v p)(x, t)f + G F V2c v n (e \x, t) , 

Here the kinematical momentum q of massless neutrino is connected with the 
canonical one, Q, as 

q = Q-G F V2c v j^(x,t) , (3) 
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and canonical definitions n — q/q — dH/dQ = dH/dq, q l = —dH/dxi = VH 
lead to (1). 

The origin of the non-conservation of the lepton current seen from (1) is the 
absence of the invariance of the Hamiltonian H (and, hence, of the neutrino 
RKE) with respect to the gauge transformation of the electron current 

j { ;\x,t)^j^(x,t)-d, X (x,t) , (4) 



where an arbitrary function x(x, t) should also obey d'Alambert equation 
d»d liX &t) = 0. 

Note that the invariance of the one-particle neutrino motion equation in a 
medium under the same gauge transformation (4) is equivalent to the neutrino 
current conservation too since in the integrand of the action S = f(L + 
L int (x, t))d 4 x there appears an additional (second) term 

Gf% v) (x, t)jj?> (f , t) - G F jfo (x, t)jj? (x, t) - G F jfo (x, t)d, X (x, t) , 

which does not contribute to the action, or such gauge transformation should 
not influence the motion equation coming from the extremum, 8S — 0, exactly 
due to djM&ij/dxp = 0. 



3 Neutrino RKE for the gauge invariant Wigner distribution 



The recipe of gauge invariance restoration is the same as in QED plasma [7]: 
we should recast (1) for the gauge invariant Wigner distribution function 

fM(q,x,t) = f^(Q,x,t) , (5) 



where the latter obeys the same RKE (1) but with the substitution of the 
kinematical momentum <f by the canonical one (3), 

df M (Q,x,t) .a/wfet) M dfM(Q,x,t) _ 

dt dx ^ dQ 



Accounting for (2), (3), the total time derivative &j( e \x,t) /&t = dj ( - e \x,t)/dt+ 
(nV)j( e \x,t) with the identity 

(nV)] {e) (x,t) - V(nJ (e) (f,t)) = -[n x V x J {e) (x,t)} , 
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and using the recipe (5) one can easily check that the RKE above takes the 
form [2] 



df^M dfM(q,x t t) (v) q» df^(q,x,t) _ 

dt +H dx + J "m^ t ) £q dq . " U > W 

where the antisymmetric tensor F^P (x, t) entering the effective Lorentz force 
is given by the weak vector current, 



F$\x,t)/G F V2c v = -V^\x,t) - t] , 

F$\x, t)/G F V2c v = e jkl {V x ]^{x , *)), , (7) 

and in accordance with (3) we changed the derivative d/dQ — > d/dq. 

Obviously, the tensor (7) (and thus the whole RKE (6)) is invariant with 
respect to the transformation (4) and obeys the continuity equation, or the 
neutrino current j^(x,t) = J(d 3 q/(27r) 3 )(q fJi /q)f^ l/ \p,x,t) is conserved, 

9£M = 0. (8) 
dx^ 



So far we did not obtain any new formulae. However, we argue that the above 
formal recipe (5) is a simple consequence of the gauge invariance under the 
transformation (4). In the next section we try to elucidate the physical sense 
of such invariance in plasma and its connection with the important definition 
of the gauge invariant Wigner distribution in SM (5). 



4 Gauge invariant Wigner distribution and induced electric charge 
of neutrino 



Let us remind the definitions and the physical sense of the gauge invariant 
Wigner distribution functions in QED plasma [8], 

f^(p,x,t) = f^(p + eA(x,t),x,t) 

= J d 3 ye^f^(x-y/2,x + y/2,t) , (9) 

where the gauge invariant distribution function in the coordinate representa- 
tion f( e \x*i, x 2 , t) is connected with the gauge non-invariant p e \xi,x 2 ,t) = 
Tr (p(t)^( e )+(f 2 )^( e )(fi)) by the important phase factor [7]: 
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f {e) (x 1 ,x 2 ,t) = exp 

x />)(£!, £ 2 ,£) . (10) 

Namely due to this phase factor the distribution (10) is invariant under the 
standard gauge transformation (with an arbitrary gauge function x(x, t) obey- 
ing the d'Alambert equation), 

tf (c) (£i) -> e- ie * (5fl>t) * (c) (£i) , 
¥ e)+ (x 2 ) -> e +iexi ' :g2 ' t) ¥ e \x2) , 
A(x 2 + £{x 1 -x 2 ),t) -> A(£ 2 + £(£i - £ 2 ),*) 

dx 2 

or, equivalently, this arbitrary phase x(zjt) cancels in (10). Such invariance 
is crucial for macroscopic physics since it provides the physical sense of the 
Wigner function (9) and the conservation of the macroscopic electric current. 

Really, as in the case of neutrino RKE (1), the kinetic equation for the gauge 
non-invariant distribution of charged particles p e \x\, x 2) t) derived from the 
quantum Liouville equation by the same Bogolyubov method does not obey 
electric current conservation. This is because the force term depends on the 
electromagnetic potentials A^x^t) which do not enter as combinations ex- 
pressed via field strengths, E, B [8]. The recasting of such RKE for the gauge- 
invariant distribution (9) allows to obtain the usual Lorentz form of the force 
term in the standard Boltzman equation for charged particles [8]: 

— dt — + v — die — e v (a: ' } [w (x ' )] J — d$ — = ' 

for which, of course, the electric current 

j( e \x,t) = J d 3 p(p fl /e p )f^ e \p,x,t)/(2Tr) 3 is conserved, dj^/dx^ = 0. 

Note that through the whole text we use exactly this standard RKE for elec- 
trons neglecting their mutual weak interactions. 

Hence, in analogy with the QED plasma definition (10) we should reformulate 
(5), which leads to the neutrino current conservation in final RKE (6), for the 
gauge invariant distribution in the coordinate representation, f^(x — y/2, x + 
y/2, t) = J dV-^/ (i %, t)/(2nf . 

To this end, comparing the neutrino canonical momentum (3) with the well- 

— * — * 

known p = P — eA(x, t) in the case of QED plasma, we find the important 



x 



(11) 
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(weak ) phase factor that connects the gauge invariant distribution f^(x±,x 2 , t) 
with the gauge non-invariant f( u \xi,x 2 ,t), 



f^\x!,x 2 ,t) =exp 



tG F V2c v {x 2 - f x ) J deJ (e) (x 2 + - f 2 ), t) 



x/W(f!,f 2 ,t) . (12) 

Our goal here is the explanation of the gauge invariance of this function via the 
gauge transformation of electron current in medium (4) with an appropriate 
electromagnetic formf actor of neutrino in plasma instead of the electron charge 
e in the standard QED transformation (11). 

Remembering the definition of the gauge non-invariant distribution f^(xi, x 2 , t) 
= Tr (^p(t)ty(^ + (x 2 )^f^ (xifj (compare with the electron case before (10)) we 
find in the completed form how the gauge invariance (4) with an arbitrary 
gauge x(x,t) is manifested in the neutrino kinetics : 



^ exp (-iGF^cvxix^t)) ¥"\x!) , 
*M+(£ 2 ) - exp (+tG F V2c vX (x 2 ,t)) *M(f 2 ) , 

•(e), 

dx(x2 + €(xi -x 2 ),t) 



r> (x 2 + ^(f 1 - x 2 ),t) -> j (e) (x 2 +C(xi-x 2 ),t) 



dx 2 

It is easy to check the cancellation of x in the phase factor in (12). 



(13) 



On the one hand, this transformation provides the gauge invariance of the 
neutrino distribution (12) and automatically the invariance of the Wigner 
function f^(q,x,t) (5) resulting in the neutrino current conservation (8). 

On the other hand, in an isotropic plasma the electric current 
jf* = — I e I if 1 i s the induced one, 

Jt\x 2 ^i{x x -x 2 \t) = J ^e^^^^M^ , 
i.e. the phase factor in (12) takes the form 



exp 



1 



GfV ~ 2CV ' 6 1 W x Jttj 0e-^«(| (£; _ l)hk 



Ana ai J V (2tt) 
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+ u,\6 ij - f ^)(e tr -l))A j (u,,k)} . (14) 

— * 

Here = (u,k) is the plasmon four-vector; ei ;tr (u>,k) are longitudinal and 
transversal permittivities in an isotropic plasma; y = x 2 — X\. 

In the two limiting cases: (i) quasistatic electric field [to <C k(v) < k, E\ — 
1 {kr D )~ 2 , u 2 (e tr — 1) — > 0) and (ii) high-frequency electromagnetic field 
(uo 3> k(v), Q 2 — * 0, e tr — 1 ks —uo 2 /uo 2 ) the phase factor (14) takes the form 
analogous to the QED plasma result (10): 



cxp 



-iel nd (x 2 - fx) j &tA(x 2 +£(x 1 -x 2 ),t) 



(15) 



where the electromagnetic field splits (see Appendix) into the longitudinal 
A{x 2 + £(xx — x 2 ),tj = A^{x 2 + £(xi — x 2 ),tj and the transversal A(x 2 + 

£(x*i—x 2 ),tj = A( tr ^x 2 +£(xi—x 2 ),tj correspondingly multiplied by: (i) either 
the quasistatic induced electric charge of neutrino or (ii) the high- frequency 
one [6], 

e md = J e|GW (16) 

_ I e | G F c v u 2 p 

2naV2 ' 1 J 



Here ro = \Jt /Anan^ is the Debye radius; u p = \J Anan^ /m e is the plasma 
frequency; T and nf 1 are the temperature and the mean electron density. 



5 The radiation damping force in isotropic medium 



Here we show that accounting for the lepton current conservation the final 
form of neutrino RKE (with an explicit dependence on the self-consistent 
electromagnetic fields E, B) does not depend whether we apply initial RKE 
(1) for the gauge non-invariant distribution function f^(q,x,t) [1], or the 
same equation written in the completed form (6) [2]. 

The situation is similar to the case of standard plasma where the initial RKE 
for the gauge non-invariant Wigner distribution /^(p, x, t) is often more suit- 
able to obtain concrete results than Boltzman equation with the Lorentz force 
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[8]. Nevertheless, the electric current should be written in the gauge invariant 
form obeying the conservation law, djjf) /dx^ = . 

Making use in (6) of the standard connection of the induced electron cur- 
rent j( e ) with the electric field E (in the Fourier representation), ji(uj,k) = 
(cu/Am) Sij{uJi k)—5ij Ej(u, k) , where the permittivity tensor in an isotropic 
plasma is given by 

e^u, k) = ei(u, k)'^pf + e tr {u, k) ( 6. 



k ^ + e tr (u,k) Uj-^r 



we easily obtain another form of the neutrino RKE (6), 



— at — + n — d~x — + e J ~l2^- t Fl{uj > k)E ^ k) 

+ F tr (cu, k)) x (E ± (u, k) + [nx B(u, k)])] ^l^A = , (1 8 ) 



where Fi(u>, k) and F tr (u, k) are the neutrino electromagnetic formfactors de- 
fined in [6], 

Fi{u, k) = Gp^cvisiiu, k) - l)Q 2 /a , 

Ftr(u, k) = G F V2c v (e tr (uj, k) - l)u 2 /a ; (19) 

Q — (u,k); E = E\\ + E±, B are the electromagnetic fields in the dispersive 
medium, and E\\ = k{kE)/k 2 . 

It is obvious that the third term on the left-hand side of Eq. (18) is proportional 
to the force of electromagnetic origin. Whereas for a point charge e (when 
the form factors are equal to unity: F\ = F tr = 1) this term is determined 
by the Lorentz force, i.e., is equal to the standard expression e(^E(x,t) + 

[nB(x, t) df(q,x,t)/dq, for neutrinos with electromagnetic structure [6], with 
allowance for the constant of the weak coupling to the electric charge, Gp ~ 
e 2 /My V , the third term in (18) is proportional to the radiation damping force 
(~e 3 ). 

The polarization origin of such a force becomes clear after simple manipula- 
tions in (18) using the explicit expressions (19) for the form factors F[ and 
F tr in an isotropic dispersive medium, for which the Fourier integrals can be 
completely calculated, and the considered term[l] 

(J- + nn ^-) PJx,t) ^M (20 ) 
e \dxjd n dtdxj J ' dqj 
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is proportional to the second derivative of the polarization vector of the dis- 
persive medium: 



which is equal to the difference between the vectors of the electric induc- 
tion, D n (x,t) = J d 4 x'e n j(x — x', t — t')Ej(x', t), and the electric field intensity 



Note that the expression (20) for the force is also valid for anisotropic media 
when the permittivity tensor depends, for example, on an external magnetic 
field. In vacuum (D = E), the effect corresponding to plasmon emission by a 
moving neutrino disappears, i.e., there is no damping force in (18). 

Finally, in special cases of the excitation in a dispersive medium of electrostatic 
waves (uj <C k(v)) or the propagation of a high-frequency transverse wave 
(a; k(v)), the term (20) can be represented in the form of the effective 
Lorentz force 



which is proportional to the neutrino induced electric charges (16) and (17), 
respectively. 



6 Conclusions 

Thus, we proved the neutrino current conservation in a medium as the con- 
sequence of the invariance of the distribution functions (12) under the gauge 
transformation (13) that is similar to the standard one (11) for the charged 
particle distribution (10). 

Second, the ponderomotive force for neutrinos given by the third term in RKE 
(6) is the damping (friction) force arising due to the neutrino electromagnetic 
structure (the electromagnetic vertex T^{uj,k) ~ GpU^ujjk)^ [6]) and is 
stipulated by the plasmon (Cerenkov) emission in medium, v — > v + 7*, for- 
bidden in vacuum [1]. 

Third, the appearance of the same neutrino electromagnetic (loop) structure 
as the factor ahead the self-consistent electromagnetic fields in the phase factor 



4nP n (x,t) = D n (x,t) - E n (x,t) , 



E n (x, t). 




dfM(q,$,t) 
dq 
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(14) looks rather natural. In particular cases of quasistatic (longitudinal) or 
high-frequency (transversal) fields there is a complete analogy with the QED 
plasma case accomplished after substitution of the electric charge e by e™ d . 
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Appendix 

In an isotropic plasma the second-quantized electromagnetic fields are addi- 
tive, = A® + A^ r \ where both the longitudinal field 

r d^k 

^ )(f ' t) = J (27r)3^(fc) (4°(fe)a (/) (fe)e~^ + g f (fe)^ )+ (A:)e^), 
and the transversal one 



obey the Lorentz gauge, dA^x, t)/dx^ = 0. Here the unit polarization vectors, 
e^e" = -1, are given by e®(k) = {k,uk)/y/^, e^\k) = (0,e< A )(A;)) , and 
for the plasmon four-vector Q M = (cu, k) obey the standard conditions for both 
components Q li e^ x \k) = 0, or fcj£^ = for the transversal part only; the 
normalization factors Nij r (k), Ni(k) = \J Q 2 <9Re £\{oj, k) / du, 

Ntr(k) = \j2uj[Etr(oj,k) + cu(9Re e~tJuj, k)/du] are given by the longitudinal 
(ei(u,k)) and the transversal (e tr (uj,k)) permittivities correspondingly. 

Note that in the Fourier representation the 3- vector parts of these four-potentials 
A^ l ' tr \uj, k) = e^ l ' tr ^ 'a(k) / 'N^ 1,1 ^ obey important opposite conditions A^(u, k) \\ 
k and kA^ tr \uj, k) = that were used above to recast Eq. (14) into Eq. (15) 
which is proportional to induced electric charges of neutrino (16), (17) corre- 
spondingly. 
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